Fingering Instability of Dislocations and Related Line Defects 
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We identify a fundamental morphological instability of mobile dislocations in crystals and related 
line defects. A positive gradient in the local driving force along the direction of defect motion 
destabilizes long-wavelength vibrational modes, producing a "fingering" pattern. The minimum 
unstable wavelength scales as the inverse square root of the force gradient. We demonstrate the 
■ instability's onset in simulations of a screw dislocation in Al (via molecular dynamics) and of a 

vortex in a 3-d XY "rotator" model. 

PACS numbers: 61.72.Lk, 62.20.Fe, 75.10.Hk, 47.32. Cc 

Modeling the motion of dislocations in crystals is mathematically challenging due to their long-iange interaction 
through. Jthe elastic strain fieldtl, and computational approaches include atomistic molecular dynamicsou and mesoscale 
modelso'El. In some contexts the complexity of the problem can be reduced by approximating the motion of an 
individual dislocation as that of a damped elastic stringtl, and such models can be used to describe the motion of 
"tin ■ other types of line defects as well, e.g. magnetic flux, lines in type-II superconductors. Recent research using this 
approach has focused on locked vs. running behavionJ of a line defect in a periodic potential due to the underlying 
. crystal lattice. A morphological instability associated with a periodic driving force at a frequency resonant with the 
periodic potential has also been studied^. 

In the present Letter we identify a fingering instability in the glide of a mobile line defect under a spatially varying 
driving force that increases along the direction of motion. This work was motivated by unpublished ma ocular 
dynamics studies of screw dislocation pair annihilation in pure Cu performed by Swaminarayan and LeSary. The 
simulations appeared to show that initially straight parallel screw dislocations of opposite Burgers vector spontaneously 
developed a waw profile prior to pair annihilation, something like the well known Crow instability of parallel vortices in 
fluid dynamicsti We were motivated by this result to consider whether the instability is driven by the attractive force 
I ' gradient between the dislocations, and indeed whether a single dislocation could be unstable under a positive gradient 
, in the driving stress. Since such stress gradients arise from common sources-e.g. external loading or dislocation 
interactions with a free surface or other microstructural features-such an instability might be active during metal 
, deformation. 

Consider the stability of an initially straight dislocation under an attractive force gradient, approximating its 
response as that of a damped elastic string. (The approximation is useful only to model the onset of any instability, 
because once the defect line deviates strongly from its initial state then nonlocal interactions can no longer be 
neglected.) The string's displacement obeys the familiar wave equation with linear damping and a spatially varying 
"p-. , driving force: 
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Here ^(z, t) is the displacement of the string in the x direction, c is the wave speed of an oscillation on the string, 7 
is a linear damping coefficient, m is the mass per unit length, and F(z,t) is the driving force per unit length along 
O ■ the string. We take periodic boundary conditions with z in the interval (0, L). For our present purposes we do not 
, include an underlying periodic potential; this is appropriate for modeling dislocation motion in FCC metals where 
? *~j ■ the Peierls stress is small and unpinned dislocations are mobile even at low temperature. 

We consider a driving force with a value that increases linearly in the direction the string is moving, as if the string 
were moving down a hill that gets continuously steeper. The force on the string at any point along its length may 
■ ' then be written: 

F(z,t)=F +F 1 V(z,t). (2) 

Any portion of the string that gets slightly ahead of the rest feels a stronger driving force. More generally Eqn. (||) 
represents the first two terms in an expansion, so that F± is defined as the local gradient of a spatially dependent 
driving force. If F± > the force gradient is destabilizing and will cause a sinusoidal perturbation of sufficiently long 
wavelength to grow rather than die away, while small wavelength perturbations are suppressed by line tension, as 
shown below. 
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To linearize the equation we switch to a (non-inertial) moving reference frame. Consider perturbations ^(z, t) from 
a stable uniform running solution ^Q(t), 

9(z,t) = 9 (t) + *(z,t) (3) 

where 

ot z m at mm 

This procedure gives rise to the linearized equation of motion for the perturbation: 

nd 2 ^ 7 <9<l Fx ~ 
Or oz z m at m 

Taking a solution of the form W = Ae ll ^ kz ~ ult \ we find the dispersion relation 



^{k) = —P-±l\^C-^(— -c 2 k 2 \ (6) 

We identify three regions of k. 

(1) Underdamped region: short wavelength (k > kd, or A < Ad) perturbations are underdamped, with 
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We identify a characteristic length Ad = 27r/fc £ z as the wavelength associated with critical damping; modes with smaller 
wavelength are underdamped. 

(2) Overdamped region: for intermediate wavelength (fc* < k < kd, or X c i < A < A*) perturbations are overdamped. 
The cut-off value k* is 



fc* = \/A- (8) 



We identify the characteristic length A* 



A* = S = W^- (9) 

k* y pi w 

Note that A* decreases as the stress gradient Pi increases, and that it is independent of the damping constant 7. 

(3) Unstable region: long wavelength modes (k < k* , or A > A*) are unstable, and will grow, leading to the 
formation of a "fingering" pattern. The fastest growing mode (which has largest positive Im u>) is always that with 
the largest unstable wavelength. 

To evaluate these characteristic lengths for a real material, we consider the case of a mobile screw dislocation in 
pure Al. The damping coefficient 7 has been estimated at about 10 -4 Pa-sedid, but in molecular dynamics simulation 
at temperatures close to zero we have found it to be about an order of magnitude smallert3. The wave speed is 
c = 29 A/psec and the effective mass depends logarithmically on system size but is in the range|-of 5 amu/A for a 
system of about 300 x 300 A 2 in cross section, again as found in a molecular dynamics simulationli 2 ]. In the absence 
of a stress gradient, the wavelength associated with critical damping at low temperature is A^ =0.2 to 0.3 /xm; at 
higher temperatures where the damping is stronger, A^ may be as small as a few hundred A. The minimum unstable 
wavelength A* depends on the local force gradient. For alumimum we estimate the energy/length of a screw dislocation 
as roughly 1 eV/A, and the driving force is a ■ b, where a is the driving stress and b is the Burgers vector=2.86 A. 
Then the characteristic unstable wavelength can be estimated roughly as A* = 5 x 10^ 4 m(a/l) 1 ^ 2 where the stress 
gradient a /I is in units of GPa/m. Thus for a screw dislocation in Al, we estimate that a stress gradient of 1 GPa/cm 
gives a minimum unstable wavelength of order A* = 50/im, and a stress gradient of 100 GPa/cm gives A* = 5/xm. 

If A* < L, where L is the length of the defect line, then all vibrational modes with wavelength A > A* will be 
unstable. While the lowest unstable mode grows the fastest, in the early stages of growth all unstable modes are 
visible and give rise to a disordered fingering structure. If we follow the evolution of the linearized equation of motion 
in Eqn. (||) beyond the early stages of growth, we observe a coarsening process in which neighboring peaks in the 
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structure merge until only a single sinusoid profile remains. In Fig. 1 we show the results of a numerical integration of 
Eqn. (H) , showing the type of patterns that arise, where the initial condition was a straight line with a small amplitude 
random perturbation. 

For dislocations in metals, the linearized equation is only valid to describe the very earliest stage of motion, so we 
would not necessarily expect to see fingering behavior like that shown in Fig. 1. However the initial growth of the 
instability should be described by the analysis above. To test this hypothesis, we performed a molecular dynamics 
simulation of the motion of a screw dislocation in a single crystal of aluminum using an Embedded Atom Method 
potential provided by Diana Farkas et alEj. To see the instability in an atomic scale simulation is not an easy task, 
because to get a minimum unstable wavelength as small as 500 A in a bulk solid, we would need an enormous stress 
gradient in the range of 10 8 GPa/m. Luckily, there is a much easier way: we make use of the force gradient produced 
by the attraction between a dislocation line and a nearby free surface. 

The force/length on a screw dislocation at a distance R from a parallel free surface scales as 1/R, so the gradient 
in the driving force diverges as 1/R 2 . When the dislocation is far enough from the surface, the local force gradient is 
small, the minimum unstable wavelength is longer than the dislocation line, and the dislocation is stable. At a critical 
distance from the wall R* (which depends on dislocation length L) the dislocation's first vibrational mode becomes 
unstable and fingers toward the free surface. As the dislocation continues to move toward the free surface, the force 
gradient gets stronger and higher order vibrational modes go unstable one by one in a cascade of instabilities. 

Our simulated system has periodic boundary conditions only along the z-axis and thus has four free surfaces; the 
dislocation is initially parallel to the z-axis which lies in the [110] direction. The system size is 45 x 45 x 183 = 370, 575 
atoms, or 111 x 105 x 524 A 3 , and the temperature is close to 0° K. The dislocation is placed a distance of 26.5 A from 
one free surface and is given a very slight sinusoidal perturbation with wavelength A = L = 524A. We estimate the 
minimum unstable wavelength as: 
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where the numerator is the energy/length of the dislocation (mc 2 ) and the denominator is the gradient of the driving 
force/lengthEl, using image terms from free surfaces ahead and behind, and canceling the energy prefactor. Here 
D = lllA is the length of the crystal in the direction of dislocation motion, and a a — b is the dislocation core size. 
For a dislocation at position R = 26. §A, Eqn. ( |l0| ) predicts A* = 258 A. Since A* is about half the size of the defect 
length L = 524 A, we expect both the first and second vibrational modes to be unstable and grow. 

The dislocation is shown in Fig. 2, where we display only atoms whose local potential energy differs from that of 
the bulk. The initial state is shown in gray, the later state in black. The dislocation splits into two partial dislocations 
separated by a ribbon of stacking fault. In each of the partials, the initial sinusoidal perturbation grows along with 
the second vibrational mode, and each of the partials fingers toward the wall. The apparent minimum unstable 
wavelength is roughly A* = 250 A, in agreement with our prediction. We have also run a simulation of the same 
set-up performed in a crystal of height 150 A; in this case the initial perturbation does not grow, as its wavelength 
A = L < A*, and the instability is suppressed by finite size. P 

Our simulation of a dislocation moving toward a free surface is closely related to Swaminarayan'sa study of disloca- 
tion pair annihilation. In our case, the dislocation is attracted to an "image" dislocation^ of opposite sign located in a 
symmetric position across the free surface. In order to see the longer time evolution of the instability, we need a larger 
system—so we have performed simulations of vortex motion in a three-dimensional XY "rotator" model of a simple 
magneto. This system is very efficient computationally, and we have run simulations of up to 4 x 10 6 spins on a 
DEC Alpha workstation. In this model, each spin on a cubic lattice has a moment of inertia 1=1 and one rotational 
degree of freedom with angular position 6i and angular velocity Wj. The spins have nearest-neighbor interactions with 
an XY Hamiltonian, an applied magnetic field B, plus a rotational kinetic energy term: 
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We derive equations of motion from Eqn. (|ll| ) under constant energy. A straight vortex line has properties analogous 
to those of a screw dislocation_in a crystal, particularly the effective mass/length, mobility behavior, Peierls barrier, 
and vortex-vortex interactionslL 2 ! 

We use periodic boundary conditions along the z direction only with four free surfaces, and size 100 x 100 x 400 
spins. In the initial state a vortex line lies parallel to the z axis at an initial temperature k^T = 0.25, in the ordered 
phase. Since the simulation is run at constant energy, fc^T fluctuates. We apply a small exteral field B=-0.003 which 
drives the vortex in the +x direction. 



3 



In Fig. 3 we show the time evolution of the vortex, starting at x=40 and moving toward the free surface at x=100. 
Using Eqn.(|o|) with a a — 1 and D = 100, we predict that the first vibrational mode should become unstable when 
the vortex reaches a distance R = 38 lattice sites from the free surface. Since this system is at finite-.temperature, the 
vortex line roughens initially due to thermal fluctuations (like those studied by Chzran and DawEj in dislocations.) 
The first vibrational mode starts to grow at about the predicted location, around x=60 (R = 40 units from the 
surface.) Around x = 65, the second mode becomes prominent (with two "peaks" visible) and grows until the vortex 
annihilates at the free surface. The second mode appears earlier than predicted by Eqn. (J10|) : perhaps nonlinear 
effects play an important role. In future work we plan to simulate a larger system in which more steps in the cascade 
of instabilities may be seen. 

The physical signifigance of this instability is an open question. Most important, we have learned how finite size 
effects play a major role in simulations of line defects in general; and that the largest characteristic length is defined by 
the smallest force gradient present. Beyond that we speculate that the characteristic length scale A* might play a role 
in the formation of dislocation patterns during work hardening. The instability could also control pairpapnihilation of 
threading dislocations in thin films, defining an effective "kill radius" that depends on film thicknessta. Confirming 
any of these mechanisms via experiment may be difficult. The motion of an individual dislocation in a thin metal 
sample can be observed via TEM, but it is not clear whether the necessary time resolution and length scales are 
technically feasible. It would be easier to observe the motion of a defect line in a liquid crystal driven by a gradient 
in an applied field, since defect motion can be viewed with an optical microscope and polarizers. 

A closely related instability has previously been identified by Z. SuoE£l who noted that a dislocation moving via 
climb can be unstable in the presence of a gradient in the density of vacancies. However since climb is a diffusive 
process, the dynamics of that instability mechanism occur on a far slower time scale where inertial effects are not 
important. However both mechanisms should play an important role in dislocation pattern formation in metals. 

We thank Lyle Levine, Robb Thompson, and Shujia Zhou for helpful discussions. This work was supported by NSF 
Grant No. DMR-59702234. 
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FIG. 1. Evolution of a line defect from numerical integration of Eqn. (10), 7 = 0.1, c = l,m = 1, and L = 1000. (a) With 
a weak force gradient Fi = 0.01, A* = 63. (b) With stronger force gradient F\ = 0.1, A* = 20. The pattern coarsens via tip 
combination. 
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FIG. 2. Molecular dynamics simulation of screw dislocation motion in Al. (a) Initial configuration, (b) Later configuration; 
dislocation has split into partials which "finger" toward the free surface. Curves are drawn to guide the eye. 
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FIG. 3. Time evolution of vortex motion in a 3-d XY "rotator" model simulation. As the vortex line approaches the free 
surface the first and then second vibrational modes become unstable. 
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